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The renormalized zero-momentum four-point coupling g r of 0(iV)-invariant scalar 
field theories in d dimensions is studied by applying the 1/N expansion and strong 
coupling analysis. 

The 0(1 1 N) correction to the /3-function and to the fixed point value g* are 
explictly computed. Strong coupling series for lattice non-linear a models are ana- 
lyzed near criticality in d = 2 and d = 3 for several values of N and the corresponding 
values of g* are extracted. 

Large- Af and strong coupling results are compared with each other, finding a good 
general agreement. For small N the strong coupling analysis in 2-d gives the best 
determination of g* to date (or comparable for N = 2, 3 with the available Monte 
Carlo estimates), and in 3-d it is consistent with available </> 4 field theory results. 
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I. INTRODUCTION. 



The study of the fixed-point behavior of quantum and statistical field theories is one 
of the central problems to be faced both from a purely theoretical point of view and for 
the purpose of investigating such phenomenologically relevant issues as the existence and 
quantitative estimate of triviality bounds. 

For understandable reasons most theoretical effort has been till now directed towards the 
analysis of a few selected models, including 0(0), 0(1), 0(2) and 0(3) in three dimensions 
and 0(A) ("Higgs") models in four dimensions. In our view, it is very useful to extend the 
analysis to the case of a generic symmetry group 0(N) and to models living in an arbitrary 
number of space dimensions d, not only for the sake of generality or in view of new possible 
physical applications, but also because this analysis may offer the possibility of testing and 
cross-checking the several different methods that have been applied to the problem at hand, 
thus putting on firmer grounds results that often rely only on a single approach and whose 
generality cannot therefore be fully understood. 

In making these statements, we have especially in mind the possibility of systematically 
extending the application of two very well known techniques: strong coupling and 1/N 
expansion. Despite their all too evident advantages, i.e. finite convergence radius for the 
strong coupling expansion and non-perturbative interpretation of results for the 1/N expan- 
sion, these techniques suffer from some drawbacks, i.e. lack of control on the accuracy of the 
resummation techniques in the strong coupling case, poor information on the convergence 
properties and technical difficulty in the extending the series for the 1/N expansion, which 
have often discouraged people from pursuing these approachs. 

Neverthless we think, and we hope to show convincingly in the present paper, that a 
renewed effort in these directions may prove very fruitful, especially when two classes of 
results are systematically compared, since agreement of 1/N expansion and strong coupling 
results for a given N can be taken as evidence of convergence of the 1/N expansion down 
to that value of iV, while the non-perturbative character of the 1/N expansion ensures the 
generality of the qualitative results (like triviality) within all the convergence domain. 

In Sec. II we introduce our notation for the general (9(iV)-invariant scalar field theory 
with quartic interaction in the continuum formulation and define the quantity which we 
are going to study, i.e. the renormalized zero-momentum four-point coupling g r , whose 
behavior in the scaling region and fixed-point value for arbitrary N and d is the object of 
our investigations. This quantity is related to the so-called Binder cumulant. We compute 
the next-to-leading 1/N correction to the renormalized coupling as a function of the bare 
coupling and the renormalized mass by reducing it to a set of Feynman integrals that can be 
evaluated in the continuum 1/N expanded model without any regularization for all d < 4. 
We show how to compute the 1/N correction to the /3-function (3(g r ) = m r and the fixed 
point value g* such that (3(g*) = 0. 

In Sec. Ill we give an exact evaluation of g r in one dimension by solving exactly the 
one- dimensional non-linear a models for arbitrary iV, and draw from this example some 
general indication about the possible dependence of g* on the parameter N. 

In Sec. IV we review the available results on non-linear a models in two and three 
dimensions on the lattice and present our explicit computations of 1/N effects for the above 
models at criticality. We briefly comment on the case d = 4 and 0(1/ N) logarithmic 
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deviations from scaling. 

In Sec. V we analyze and discuss the strong coupling series for the renormalized coupling 
of the non-linear a models in two and three dimensions on the lattice and for arbitrary iV, 
which we extracted from the results of Ltischer and Weisz [1] as elaborated by Butera et 
al. [2]. Strong coupling results are compared to all available calculations presented in the 
literature (</> 4 field theory at fixed dimensions and Monte Carlo simulations) and to our 1/N 
results, finding a good general agreement. 

Finally in Sec. VI we draw some conclusions. 

An explicit representation of the 1/N correction to the /3-function is exhibited and di- 
scussed in the Appendix. 



II. THE RENORMALIZED COUPLING AND ITS 1/N EXPANSION. 

According to the previous discussion it is interesting to form a renormalization-group 
invariant dimensionless combination of vacuum expectation values playing the role of a 
renormalized four-point coupling and to study its behavior in the proximity of a critical 
point. In particular we are interested in (9(A^)-invariant scalar field theories in arbitrary 
dimensions d < 4 and we wish to apply 1/N expansion techniques to the above-mentioned 
problem. 

From the point of view of the 1/N expansion the standard notation is somewhat incon- 
venient: we shall therefore define our own conventions, trying to establish correspondence 
with the literature as far as possible, and especially trying to make all relationships with 
Refs. [1,3-6] as transparent as we could. 

The usual (9(A^)-invariant Euclidean continuum Lagrangian takes the form 

c = ^JdJ+^ + ^y. (i) 

It is however convenient to redefine the quartic coupling (both bare and renormalized) 
according to the definition 

Ng 

(2) 



We shall also define 
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go 

l 

7 = — , 

go 



(3) 



and introduce an auxiliary field a. 

The resulting effective Lagrangian is 



N r i 

— [Pd^s + i[3a{s 2 -l) + 7 a 2 ] , (4) 
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and after performing a Gaussian integration over the field J* we obtain 



N r ^ 
C = — [Tiln (3 (-d^ + ia) - i(3a + 7a 2 J , 



(5) 



which reduces to the usual effective large- N action for the non-linear a model in the limit 
7^0. 

Correspondence with Refs. [1,3-6] is established by the relationships 
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2k, 
iV ' 

N2X 



(6) 



Renormalization is performed according to the following prescriptions for the two and 
four-point correlation functions of the field <f>: 



r (2) (p,-p) a/3 = r( 2 )(p 2 )^, 

r (2 V)= z;' [ml + p* + o(p*)\ , 



(7) 



anc 



rW(o) : 

Let's now notice that 



7-2 9r / 2\2- 



,r (4) (o,o,o,o) c 



-N ~ V'-'-'-Wr ( 2 )2 -f - 

[r( 2 )(o,o) aa ] 2 K r) V n 



i + 



(8) 



(9) 



Equation (9) will be our working definition of the renormalized four-point coupling. 

In order to compute the leading and next-to-leading contributions to g r in the continuum 
1/N expansion, we shall need an evaluation of the corresponding contribution to the two- 
point function and to the zero-momentum four-point function. 

The evaluation of the Feynman rules shown in Fig. 1 is essentially straightforward. We 
only mention that the bare propagator of the <f> field is expressed in terms of a "bare" large- N 
mass parameter ra 2 , introduced by the gap equation 



/ 



d d p 



(2tt)<V + m 2 



P + 2 7 ™ 2 



(10) 



while the propagator of the Lagrange multiplier field a is defined to be -^A^ra 2 ,), where 
in turn 



d d p 1 1 
(2ir) d p 2 + ml (p + k) 2 + ml 



+ 7 



±£(2-1) /*» J 

2 (4tt)^/2 \A + % 



o _ i I i 

2'2'2'l + k 2 



+ 7 



(11) 
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The relevant higher-order lagrangian effective vertices are obtained by taking derivatives of 
A -1 (A;) with respect to ra 2 ,, according to the correspondence table 



VW(0,k,k) 



d 



dml 



A~\k) 



2V^(0,0,k,k) + V^(0,k,0,k) = _-£-V(3)(o,fc,fc) = ^—A" 1 ^), 



dml 



(dml) 



(12) 



where mass dependence is suppressed in the arguments. The derivatives appearing in 
Eq. (12) may be evaluated by a generalization of the so-called "cutting rule" of Ref. [7], 
whose d- dimensional form is 



d 



dml 



A~\k) 



k 2 + Ami 



[(3 - d)A-\k) + A-^O) + 7 (<Z - 4)] 



(13) 



In writing Eqs. (10) and (11) some ultraviolet regularization, when needed, is assumed. 
Actually our final results will turn out to be independent of the regularization as expected 
on physical grounds. 

Equation (9) shows that in order to compute g r to any definite order in the 1/N expansion 
we must be able to compute the quantities T^ 2 \p), r( 4 )(0) and m 2 with the same precision. 
Leading order calculations are straightforward. Next-to-leading contributions may be for- 
mally represented in terms of a few fundamental integrals, which are graphically represented 
in Fig. 2 and listed below: 



d d k A(k,m 2 ) 



(2ir) d (p + k) 2 + ml 



2 ' 



(14) 



4 6) K) = -U(0,ml) I -0^V^(O,k,k)A(k,ml) 



(15) 



B ^ = h 



d d k A(k,m 2 ) 2 



(16) 



(2iv) d (k 2 + ml) 2 ' 

It is very easy to show that the two-point function and the renormalized mass are respectively 



r( 2 V) = p 2 + ml + 1 [s(*y X) + + o (^) 



(17) 



anc 



m l = m o + 



N 



sW(0,mS) + Si°VS)-mS 



dp 2 



(a) 



(o,m 2 ) 



+ 



(18) 



Explicit use of Eqs. (12), leading to the graphical identities drawn in Fig. 3 allows to obtain 
the representation 
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-ivrW(o) = A(o , m n 



i + 4A4 a) (o,^ ) + l-^4 6) K) 



N dm 2 



N dm 2 ' 



2 

N 



^ Bl (ml) + 0(—^ 



1 

\N 2 J 



(19) 



It is now important to notice that, in order to obtain a finite result, g r must be expressed 
in terms of the renormalized mass m 2 . This is achieved by inverting Eq. (18), which leads 
to 



2 2 

m n = m. 

r N 



Si a) (0,m r 2 ) + S^K)-m 



dp 2 



(0,m r 2 ) 



(20) 



and as a consequence 



r( 2 V) 

and 



sW(p 2 ,m r 2 )-Sr(0,m r 2 ) + 



(»)/ 



m; 



,gsj a) 

<9p 2 



(o,^ 2 ) 



+ 



(iV 2 ) 



(21) 



A(0,m 2 ) = A(0,m 2 ) + ^A(0,m 2 )[l- 7 A(0,m 2 )] ^ - 2^ 



S^(0,m 2 ) sj b V 2 ) as<"> , x 



m: 



m; 



<9p 2 



(22) 



Collecting all the above results and substituting into Eq. (9) we obtain the following 
representation of the renormalized coupling: 



(m 2 )f- 2 A(0){l+l[l- 7 A(0)] (^-2 



+ 



mi 



mi 



dp 2 



(0) 



+ 



N 



2 ^M + ^_ 2 ^ (0) _ 2 A-(0) Bl 



dm 2 



dm 2 



dp 2 



(23) 



where all quantities are now computed with m 2 , replaced by m 2 and mass dependence is 
suppressed in the arguments. Substituting Eqs. (14), (15) and (16) and making explicit use 
of Eq. (13), one obtains the following representation: 



g r = (m 2 )- 2 A(0) + - [1 - 7 A(0)] (3 - d)2' 



d-i 



3m 2 



( f-4 2(^-1) ' 

(k 2 + m 2 ) 3 ^ (A; 2 + m 2 )(k 2 + 4m 2 ) \k 2 + m 2 k 2 + 4m 2 



2m 2 



1 r ddk a/,x 

1 / , ^ a /" ^ a 2(<Z — 1) 



+ 



3m 



2(d — 1) 

(F + m 2 ) 3 ' (F + m 2 )(k 2 + Am 2 ) \k 2 + m 2 ' A; 2 + 4m 2 



+ 



2tt)° 
<Z d fc A(A;) 2 
27r) d (A; 2 + 4m 2 ) 2 



(A; 2 + m 2 ) 



(d-4) 7 A(0) 



+ [W 2 



(24) 
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Further computational simplification is achieved by making use of the following straightfo- 
rward consequence of Eq. (13): 

A-mMk? = *m; + (4-<Q* 2 A(h) _ 2P(P + 4^) dA(k) 

K > K > 4m 2 + (4 - d) 7 A(0)P ^ ; 4m 2 + (4 - d) 7 A(0)P dk 2 ' ^ ; 

which may be applied to Eq. (24) in order to get rid of the A(&) 2 dependence in the integrand, 
while a partial integration may eliminate the dependence on d gffi • It is easy to recognize 
that whenever d < 4 and 7 > all integrations are finite. The final result can be formally 
expressed by the relationship 

9r = 9l°\x) + ±g?\x) + (^) , (26) 

where all dependence on the renormalized mass and the bare coupling can only come through 
the dimensionless combination x = m^~ d /g . Specifically one obtains 

5S 0) = ( 1 + *) 1 , (27) 



where 



1 m 



4-<f 



g. 2 J (2ir) d (p* + m 2 ) 2 2 (47r) d / 2 ^ ; 

is the inverse of the large- N fixed point value of the renormalized coupling. 

Equation (26) is the obvious starting point for the construction of the 1/N expanded 
/3-function of the model, via the relationship 

m) = m rjt r = (4 - ^fr'-(-) = + iv^ 1 ^ + ' (29) 

where x(g r ) is obtained by inverting the equation g r = g r (z), and it admits in turn a 1/N 
expansion in the form 



x 



(g r ) = x Wfa) + ± x Wfa) + o(jp) . (30) 



It is easy to recognize that Eq. (27) implies 



X {0 \dr) = ^ " - , (31) 
9r 9* 

and therefore the large- N limit of the /3-function reduces to 

^°\g r ) = (d-4)g r [l-^j . (32) 
This is the standard (large- N) one-loop result provided that we identify 
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A-d 



r( 3 -f) 



g* (47r) d / 2 

A simple consequence of Eqs. (26), (27), (30) and (31) is the relationship 



(33) 



(34) 



We may now consider the expansion of Eq. (29) in powers of 1/N and notice that the 
derivative of g^\x) with respect to x^ may be exchanged with a derivative with respect to 
g r . As a consequence after some manipulations we can prove the relationship 



^\g r ) = [/?«»(&)] 



2 d 



gi 1] (* (0) (g,)) 



(35) 



where notable simplifications occur when evaluating g^ directly as a function of g r 1 — g t 1 . 
In particular 



m„ 



d-4 



A(k) 



1 + g r IL(k/m r ) 



where 



= 2/(2? 



d d p m 



4-<f 



(2tt)<V + 



mi 



(p + k) 2 + ra 2 p 2 + ml 



(36) 



(37) 



is a regular dimensionless function with the property 11(0) = and a finite d — > 4 limit. 
As long as d < 4 one may show that ^ x \g r ) is well-defined and finite for all < g r < g*. 
We obtained an explicit integral representation of ^ x \g r ) for arbitrary d, and showed that 
the series expansion of such a representation in the powers of g r may be obtained also in 
the d — > 4 limit and reproduces all known results as long as comparison is allowed. The 
representation of ^ x \g r ) and a short discussion of its features are presented in the Appendix. 
The non-perturbative properties of f3^\g r ) when cZ — > 4 will be analyzed and discussed in a 
separate publication. 

For what concerns the very important issue of the fixed point of g r , we must notice that 
the /3-function vanishes when x — > 0, i.e. when 

g* r =g* + -^(o) + o (^) . (38) 

Equation (24), supplemented with Eq. (25), lends itself to an easy evaluation in the limit 



x — > 0, corresponding to the limit 7 — > 0. The final result is 

d d k Ao(Jfe) / 4m 2 



<^(0) 
where 



( 3 ^ 2 / (2ir) d (k 2 + m 2 ) 2 \k 2 + ™ 2 



d 



m~ 



k 2 + 4ra 2 



(39) 
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A (k) = limA(Jfe) 

x— »0 



(40) 



Notice that the fixed-point value of the renormalized coupling may be obtained directly 
by computing the 7 — s- limit of the coupling g r in the scaling region. However this is 
nothing but the value taken by g r in the corresponding continuum limit field theory, that 
is the usual non-linear a model in d- dimensions. In turn this is the limit of the lattice 
non-linear a model when /3 — > /3 C , the value of the coupling such that the renormalized mass 
(i.e. inverse correlation length) is equal to zero in the lattice 7 — s- limit. 



III. NON-LINEAR a MODELS IN ONE DIMENSION. 

Before discussing the general d- dimensional case, let's illustrate some features of the pro- 
blem by solving the simple but not trivial one dimensional case. One-dimensional non-linear 
a models are a completely integrable system, both on the lattice and in the continuum [8-10]. 
Indeed in any lattice formulation with nearest-neighbor interactions the two and four-point 
correlation functions are easily expressed in terms of two quantities that in turn are related 
to vacuum expectation values of the model defined on a single link. 

Without belaboring on the rather trivial manipulations needed to derive these re- 
sults [8,9], we simply quote that, in any (9(iV)-invariant a model theory satisfying the 
constraint s-s = 1, one may write 



( On ) 



T>n—m cab 



S m S m 



for n > m, and 



cab 



N 



S n S n 



ccd 



N 



rjn—m 
B 22 



cabcd 



ab ccd 



N(N + 2) 



N 2 



-8 ab & 



(41) 



\ rr 



„b c d \ 
S„S„S„ ) c 



m n p q 



n—m-\-q—p 



cabcd 



Tjp-n 
#22 



R 2(p-n) N 
n ll 



N(N + 2) N 2 



(42) 



for q > p > n > m, where m,n,p, q are integer numbers labelling lattice sites, 

^abcd ^ab^cd _j_ ^ac ^bd _j_ ^ad^bc 

and 



(43) 



£11 
B22 



(si - s ) , 
N( (i*! • s ) 2 ) 
N - 1 



(44) 



where expectation values are taken in the single-link model: Bn and B 2 2 are the character 
coefficients in the (pseudo) character expansion of the model [8,9], or, equivalently, the coef- 
ficient of the expansion of the theory in hyperspherical harmonics. The results corresponding 
to different orderings of the lattice points are obtained by trivial permutations. 
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Zero-momentum lattice Fourier transforms can be computed as functions of Bu and B 22 
by performing trivial summations of geometric series. One may then easily recognize that 

T (2) (0) aa = (45) 
1 — #11 

and 

, (1 - B u ) 2 , , 

m\ = 1 11} , (46) 

-Dll 

while after some purely algebraic effort one obtains 

As a consequence by applying Eq. (9) one obtains 
2 



1 + 



2(1 -Buf ANBt, 16gn n l-^n „ ' 



(48) 



We now want to take the critical limit, which, as shown in Ref. [10], depends in a complicated 
way on the specific hamiltonian, as in one dimension there are infinitely many universality 
classes. We will restrict ourselves to those theories for which m r — > 0. Thus Eq. (46) implies 
Bu -> 1. When N > 1 we have also [10] B 22 ^ 1, so that 

l + jf)X =2(N + 2)-A(N- 1)b ^\^-. (49) 

The quantity \ZbII characterizes the universality class being simply the ratio of the mass- 
gap in the spin-one and spin-two channels. Within the universality class corresponding to 
the standard continuum limit, we have 

1 - Bu N - 1 

lim TT 1 = 7T- ( 50 ) 

Bu-i 1 - B 22 2N K ' 



and as a consequence 



1 + = *{ 1 -ir) (51) 



for N > 1. This solution agrees perfectly with the prediction resulting from Eq. (39). When 
N < 1 we have no general argument for the behavior of B 22 in the massless limit. We 
may however restrict our attention to the universality class corresponding to the standard 
continuum limit, and within this class we may consider the specific lattice example of the 
minimal nearest-neighbor coupling. For this action one may show that for arbitrary N 

I N/2 (NB) 

-Dll ~~ 



I N /2-l(N3) ' 



d _ In,2+i(NB) _ Bu 

B22 ~ iw-W) - l ~T' (52) 
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When N > 1 B\\ is strictly smaller than one for all finite values of (3 and only in the limit 
(3 — > oo the massless regime is attained, in which case B 22 -> 1 as well, as expected from 
the general argument. However when N < 1 one may numerically check that a finite value 
(3 C exists such that Bn((3 c ) = 1: as a consequence B 22 ((3 C ) = 1 — 7^ 1. Since B22 is strictly 
different from one in the massless limit, we get from Eq. (49) that within this universality 
class 

(l + ^) T T = 2(N + 2) (53) 

for N < 1. Let's notice that the two solutions connect very smoothly to each other (the 
function and its first derivative have the same left and right limit) because of the double zero 
at N = 1 in the contribution that does not vanish when N > 1. Fig. 4 shows (l + g* 
versus N. 

This is very reminiscent of what is going to happen when d = 2: for small N f3 c is finite 
and there is a domain of analiticity in N for g r around N = —2, while for large N f3 c is 
infinite (asymptotic freedom) and analiticity in 1/N is present. The two regimes seem to 
meet smoothly at N = 2. 

As a further check of our results we may consider the N = case. This is a very simple 
model of self-avoiding walks in one dimension. All computations are straightforward and 
one obtains, with the notations adopted here, 

fc( 1 + f)s = 4 . < 54 > 

in full agreement with our general formula. Also intermediate steps are reproduced, with 
the identification B\\ = (3, B 2 2 = 0. 



IV. NON-LINEAR a MODELS IN HIGHER DIMENSIONS. 

Lattice non-linear a models, which we may choose to describe in terms of the standard 
(9(iV)-invariant nearest-neighbor action 

Sl = -Nj3^2s x ■ s x+ll (55) 

subject to the constraint i*J = 1, when considered on a d- dimensional lattice with d < 4 have 
a nontrivial critical point (3 C < 00 whose neighborhood (scaling region) is properly described 
by a renormalized continuum field theory. This theory is in turn nothing but the 7 — )• limit 
of the standard (9(iV)-invariant scalar field theory (linear a model). We may therefore study 
the critical properties (and in particular the fixed-point value of the renormalized coupling) 
of the symmetric phase of the O(N) model by exploring the region /3 — > /3 C of the lattice 
model. 

The left-hand-side of Eq. (9) has a simple reinterpretation in terms quantities defined 
within the associated lattice spin model. Setting 
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X 

m 2 = ^x 2 ( s ■ s x ) , 

X 

£2 m 2 1 



2^X ' 

%4 = XI ( ^0 • Sx S y ■ S z } c , (56) 



one can argue that the combination 

X4 



x 2 i d 



(57) 



should either admit a non-trivial limiting value or vanish with logarithmic deviations from 
scaling when the critical line is approached. This is essentially a consequence of the existence 
of an unique diverging relevant scale in the scaling region. It is furthermore trivial to show 
that in the scaling region, m r — > 0, 



and in particular 



/* = /(&) = + (59) 



We also mention that / can be written in terms of the Binder cumulant defined on a L d 
lattice 



2 ({S-S) 2 ) 



U L = 1 + \ J 2 , (60) 

N {(s-s)) 



where S = J2 X Indeed 



f = "HHL^iff ■ < 61 » 

As already mentioned in Sec. Ill, there is a crucial dependence on the space dimensio- 
nality as well as on N. In two dimensions it is well known that models with —2 < N < 2 
are well described at criticality by conformal field theories with c < 1. In particular at 
N = —2 the fixed point is Gaussian, N = corresponds to a model of self-avoiding ran- 
dom walks, N = 1 is the solvable Ising model, and N = 2 is the XY model showing the 
Kosterlitz-Thouless critical phenomenon, characterized by an exponential singularity at a 
finite (3 C . 

When — 2 < N < 2 the critical point occurs at a finite value of (3 C , which should coincide 
with the convergence radius of the strong coupling series. When N > 3 there is apparently no 
criticality for any finite value of (3. This is consistent with the Mermin- Wagner's theorem on 
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the absence of spontaneous symmetry breakdown for two-dimensional continuous symmetry 
and with the weak coupling (large (3) prediction of asymptotic freedom and dynamical mass 
generation for this class of models. Large- N results and the 1/N expansion are completely 
consistent with the above picture. From the point of view of the renormalized coupling 
analysis it is however impossible to distinguish between the two behaviors, since they are 
both compatible with a non-zero value of /*. 

We now briefly present some large- N results regarding / and its limit at (3 C . On the 
lattice, using the action (55), the large- N limit of f(f3) can be easily obtained from the 
saddle-point equation 



I 



d d q 



2 ' 



(62) 



where q 2 = 4 sin 2 (g^/2), and the relation 

m = ■ 



4-d a/3 

o dm 2 



(63) 



In 2-d the above equations are made more explicit by writing 



p = —kK(k) , 



(64) 



(2 \ —1 
1 + ^) , K and E are elliptic functions. Fig. 5 shows f(f3) versus (3. In the 

large-/? limit the continuum result (28), i.e. /* = 8"7r, is recovered. 

Our 1/N expansion analysis of Sec. II leads to the evaluation of the 0(1 /N) correction 

to /*, indeed in two dimensions: 



r 



8tt 



1 + l± + o(— 

N \N 2 



(65) 



with j x = -0.602033.... 

In three dimensions we face a quite different situation. The critical point occurs at a 
finite value (3 C for all values of N. At N = oo [11] 



(3 C {N = oo) 



t d 3 q 1_ 
J (2tt) 3 J 2 



0.252731. 



(66) 



The four-point renormalized coupling f(f3) at N = oo is shown in Fig. 6. At the critical 
point /* = 16-7T. It is also possible to estimate the deviation of (3 C (N) from Eq. (66) by the 
1/N expansion technique presented first in Ref. [11], leading to the relationship 



(3 C (N) = &(oo) + | + O (^) 



(67) 



where 
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(68) 



Numerically h\ ~ —0.117. It is also important to possess an estimate of the value of the 
internal energy E at the criticality 



E C (N) = E c (oo) + ^ + (±) 



(69) 



We notice that 
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(70) 



where S^ra 2 ,) is the lattice counterpart of S^ra 2 ,), defined in Eq. (15). By setting d = 3 
and by considering the ra 2 — > limit we then obtain 



E C {N) 



+ 



6/3 c (oo) 6iV 



6x 2 

+ 



&(oo) 2 &(oo) 7 (2tt) 3 q 2 



y ?27r 
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(71) 



Numerically we obtained E c (oo) = 0.340537... and e\ ~ —0.07. 

For the three dimensional case our 0(1/ N) calculation of /* gives 



/* 



16?r 



N \N 2 J 



(72) 



with /x = -1.54601.... 

In the two and three dimensional O(N) models a number of techniques have been applied 
to the determination of /*. In particular we mention the </> 4 field theoretical approach at 
fixed dimensions proposed by Parisi [12] and developed in Refs. [13,14], making also use of 
Borel resummation techniques (see for example Refs. [15,16] for a review on this approach). 
This method has been applied to N = 1 (Ising models) in 2-d [14], N = 0,1,2,3 [13,14,17] 
and many larger values of N [17] in 3-d, leading to rather precise estimate of /*, especially 
in 3-d. In order to compare our results with the field theoretical calculations we must keep 
in mind that it is customary to rescale the coupling in such a way that for all values of N 
the one-loop fixed point value of the new coupling g be exactly one [13,14]. By comparing 
the one loop expression of the /3-function, which in our notation would be 



P(g r ) = (d-4:)g r + 



N + 8 T{3-d/2) 
_ N (47r) d / 2 



t + 0(3?) 



(73) 



we 



find 
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iV + 8r(2-^/2) 
^ iV + 2 2(47r) d / 2 y K ' 

We also mention a determination of /* for N = in 3-d by working directly with the 
self-avoiding random walk model [18], which turns out to be in full agreement with the 
corresponding </> 4 field theoretical calculation [13,14]. 

Estimates of /* can also be obtained by Monte Carlo simulations using the lattice for- 
mulation of the theory, by directly measuring /(/?)• Numerical studies concerning the four- 
point coupling have been presented in the literature for some two-dimensional models: for 
N = 1 [19], and N = 2,3 [20]. The comparison with these works must take into account the 
extra factor N in our definition (61) of the four-point coupling /. 

Finally let's briefly comment on the d = 4 case. In this case it is not possible to define 
a non-trivial limit for the non-linear a model in the strict 7 = regime, at least within the 
1/N expansion, since we obtain the naive result A(&) = implying g* = 0. This is however 
consistent with the common expectation that (9(iV)-invariant models in four dimensions 
may only have a trivial fixed point, in which case the critical region should be characterized 
by logarithmic deviations from scaling. That this is the case has been shown by Liischer and 
Weisz by making use of the strong coupling expansion in a beautiful series of papers [1,3-6]. 
Kristjansen and Flyvbjerg [21] in turn have developed the lattice 1/N expansion of the O(N) 
invariant models in four dimensions both in the symmetric and in the broken phase, finding 
substantial agreement with Refs. [1,3-6] at N = 4 in the region around criticality. 

We may add that, by properly manipulating the expression presented in the Appendix 
in the limit c? — > 4, it is possible to compute exactly (albeit only numerically) the 1/N 
correction to the /3-function of the (</> 2 ) 2 model for all values of the (running) renormalized 
coupling. In turn these results might be used to improve our understanding of the non- 
perturbative limit of a strongly interacting Higgs sector on the line traced by Refs. [22,23], 
where the leading order result was analyzed. 



V. STRONG COUPLING ANALYSIS. 

The non-triviality issue can be also investigated by high-temperature series methods 
formulating the theory on the lattice. We consider the nearest-neighbor formulation (55) of 
O(N) vector models. Notice that in (55) we have introduced a rescaled inverse temperature 

The strong coupling expansion of f(/3) has the following form 

(75) 

Series up to 14 th order of the quantities involved in the definition of /(/?), i.e. %, ra 2 and 
% 4 , have been calculated by Liischer and Weisz [1], and rielaborated by Butera et al. [2]. 
From such series one can obtain A d ((3) = (3 d / 2 f X/3) up to 13 th order. We mention that, 
for other purposes, we have calculated the strong coupling series of the Green function 
G(x) = ( s ■ s x ) up to 21 th order in 2-d, and up to 15 th in 3-d, obtaining the corresponding 
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series of the energy, the magnetic susceptibility and the second moment correlation length 
to the same order [24]. 

We also considered strong coupling series in the energy f(E), which can be obtained 
by inverting the strong coupling series of the energy E = (3 + 0(f3 3 ) and substituting in 
Eq. (75): 

(76) 

From our strong coupling series of the energy [24], we could calculate B d (E) = E d l 2 f(E) 
up to 13 th order. 

Before describing our analysis of the above series based on Pade approximants (PA's) 
technique (see Ref. [25] for a review on the analysis of strong coupling series), we recall that 
PA's are expected to converge well to meromorphic analytic functions. More flexibility is 
achieved by applying the PA analysis to the logarithmic derivative of the strong coupling 
series considered (Dlog-PA analysis), and therefore enlarging the class of functions which 
can be reproduced to those having branch-point singularities. In general more complicated 
structures may arise, such as confluent singularities, which are sources of systematic errors 
for a PA (or Dlog-PA) analysis. In particular confluent singularities at /3 C , i.e. confluent 
corrections to scaling arising from irrelevant operators [14,26], lead in general to a non- 
diverging singularity of f([3) at f3 c . Indeed in the presence of confluent singularities we 
would expect f([3) to behave as 

f{P) ~ /* + c(&-/3) A (77) 

close to (3 C with A > 0. Such a behavior close to (3 C cannot be reproduced by a PA's or 
Dlog-PA's, while it could be detected by a first or higher order differential approximant 
analysis [27]. Therefore in order to reduce systematic errors one should turn to more general 
and flexible analysis, such as differential approximants, which, on the other hand, require 
many terms of the series to give stable results. We tried also this type of analysis without 
getting stable and therefore acceptable results, very likely due to the relative shortness of the 
available series. We then expect PA and Dlog-PA analysis to be subject to larger systematic 
errors when confluent singularities are more relevant, as in 3-d models at small iV, where 
they represent a serious problem also in the determination of the critical exponents from the 
available strong coupling series. 

It is important to notice that the accuracy and the convergence of the PA estimates may 
change when considering different representations of the same quantity, according to how well 
the function at hand can be reproduced by a meromorphic analytic function. By comparing 
the results from different series representations of the same quantity one may check for 
possible systematic errors in the resummation procedure employed. To this purpose, in our 
study we will compare estimates of /* coming from the strong coupling series of both f(/3) 
and f(E). Exact results at N = oo presented in the previous section, beside giving an 
idea of the behavior of f(/3) at finite iV, represent useful benchmarks for strong coupling 
methods. 

The most direct way to evaluate /* = f(f3 c ) would consist in computing [l/m] PA's 
■Ai/m{fi) from the available series of Ad(/3), and evaluating /3~ d / 2 yl;/ m (/3) at the critical point 
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(3 C (at least if (3 C < oo; if (3 C = oo things are trickier as we will discuss below). This simple 
procedure works already reasonably well, but we found more effective a Dlog-PA analysis, 
which showed a greater stability and whose results will be presented in the following. 

Our Dlog-PA analysis consisted in computing [l/m] PA to the strong coupling series of 
the logarithmic derivative of Ad(f3), let's indicate them with ~D\og^ m Ad(/3), and then a set 
of corresponding approximants fi/ m (/3) to /(/?), which are obtained by reconstructing f(f3) 
from the logarithmic derivative of Ad((3): 

hUP) = ^/i ex P J*dP'T>log l/m MP') , (78) 
All the approximants with 

Z + m > 10 , m > Z > 4 . (79) 

were considered in order to check the stability of the procedure. Notice that, for given Z,ra, 
the number of terms of the series of Ad((3) used by the corresponding Dlog-PA is n = Z+ra+1. 
Once the approximant fi/ m (/3) is computed, if (3 C is finite, its value at (3 C provides an estimate 
of /*. 

This requires a rather precise determination of (3 C , which is in some cases available in the 
literature from strong coupling and numerical Monte Carlo studies. When (3 C was not known, 
we estimated it from a Dlog-Pade analysis of the strong coupling series of the magnetic 
susceptibility (up to 21 th order in 2-d and 15 th in 3-d). Our strong coupling determinations 
of (3 C in 3-d models at large N compare very well with the 0(1 /N) calculation (67), as shown 
in Fig. 7. Let's notice that the error on the value of (3 is small enough not to be relevant for 
the estimate of /*. 

In order to better understand the analytic structure of f(f3) we have done a detailed 
study of the complex-plane singularities of the Dlog-PA's of Ad(f3). We have first checked 
hyperscaling. A violation of hyperscaling would lead to a behavior f(f3) ~ Ad(f3) ~ (f3 c — f3) e 
for (3 — > (3 C , and thus the Dlog-PA's would show a simple pole at (3 = (3 C . We recall that a 
Dlog-PA analysis is in general very efficient in detecting power-law singularities. We have 
found no evidence of such a pole, confirming hyperscaling arguments. However notice that 
Eq. (77) when 1 > A > implies a behavior 

DlogA d (/3) ~ (/3 C -/3) A - 1 (80) 

close to (3 C . In two dimensions A ~ 1 and therefore we do not expect to find singularities 
around (3 C . This is confirmed by the analysis of A 2 (f3). In three dimensions instead, at 
least for small iV, A ~ 0.5 is expected, thus A 3 (f3) should behave as in Eq. (80), and in 
the Dlog-PA's the singularity should be mimicked by a shifted pole at a (3 larger than (3 C . 
Indeed in the analysis of A 3 (f3) we have found a singularity typically at (3 ~ 1.1 -f- 1.2 (3 C . 
This fact will eventually affect the determination of f(f3) close to (3 C by a systematic error. 
However since the singularity is integrable the error must be finite, and the analysis shows 
that such errors are actually reasonably small. 

Sometime PA's showed spurious singularities on the positive real axis (or very close 
to it) for (3 ^ (3 C . We considered these approximants defective, and discarded them from 
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the analysis. Such defective PA's were a minority, as the Tables show. The only stable 
singularity detected by the Dlog-PA's of Ad((3) lies in the negative (3 axis and closer to the 
origin than (3 C : it turns out to be nothing but a regular zero of Ad((3). The position of this 
negative zero is reported in the Tables I,II,III and V for several values of N. 

As final estimates of /*, reported in the Tables IV and VI, we take the average of the 
values fi/m((3c) from the non-defective PA's using all available terms of the series, i.e. those 
with n = Z+ra+1 = 13. The errors displayed in the Tables IV and VI are just indicative, they 
are the variance around the estimate of /* of the results coming from all PA's considered 
(cfr. (79)), which should give an idea of the spread of the results coming from different 
PA's. Such errors do not always provide a reliable estimate of the systematic errors, which 
may be underestimated especially when the structure of the function (or of its logarithmic 
derivative) is not well approximated by a meromorphic analytic function. In such cases a 
more reliable estimate of the systematic error would come from the comparison of results 
from the analysis of different series representing the same quantity, which in general are not 
expected to have the same structure. 

For this reason we have considered the series in the energy variable, which we have 
analyzed exactly as f{(3). In this case instead of (3 C we needed E c , the energy at the critical 
point. When the value of E c was not available in the literature, we estimated it by the first 
real positive singularity found in the analysis of the available strong coupling series of the 
magnetic susceptibility expressed in powers of E. This procedure provides an estimate of 
E c much less precise than (3 C (see Table V for the values obtained in 3-d), but sufficiently 
good to our purposes, given that f(E) is smooth around E c . In Fig. 8 we compare our 
determinations of E c in 3-d models with the large- N result (71), showing agreement within 
the uncertainty of the strong coupling results. 

In asymptotically free models where (3 C = oo, the task of determining /* from a strong 
coupling approach appears much harder. On the other hand, since at sufficiently large (3 we 
expect that 

W)- r I ~ r 2 , (si) 

a reasonable estimate of /* could be obtained at /3-values corresponding to large correlation 
length, £ > 10 say, where the curve f((3) should be already stable (scaling region). Notice 
that this is the same idea underlying numerical Monte Carlo studies. Another interesting 
possibility is to change variable from (3 to the energy E, and analyze the series in powers 
of E. In the energy variable the continuum limit is reached for E — > 1, and therefore the 
strong coupling approach to the continuum limit appears more feasible. In order to reach 
the continuum limit from strong coupling, we believe this change of variable to be effective 
especially for the analysis of dimensionless ratios of physical quantities. 

To begin with we present the results obtained for the 2-d models. PA's for — 1 < N < 2 
are quite stable, giving estimates of /* very close to each other as shown in Table I, where the 
values of the approximants fi/ m (/3) at (3 C are reported. For N < —1, due to the instability 
of the corresponding PA's, we could not get reliable estimates of /*. Fig. 9 shows fe/6((3) 
versus (3 / (3 C for various values of N . Differences in the other PA's were of the order of the 
width of the lines drawn in the figure. Final estimates of /* are reported in Table IV. In 
order to check possible systematic errors in our analysis, we applied the above procedure to 
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the 13 th order series of B 2 {E) of the Ising model (N = 1). The value of f(E) at E c = ^2/2, 
the energy value at (3 C , must give again /*. As in the analysis of A 2 ((3), approximants 
fi/ m (E) turn out to be rather stable, leading to the result /* = 14.64(10), which is perfectly 
consistent with the estimate /* = 14.63(7) coming from the analysis of A 2 ((3). 

For the Ising model (N = 1) our estimate of /* is in agreement with the result of 
Ref. [28], obtained by a slightly different strong coupling analysis (the value reported there is 
/* = 14.67(5)), and with the estimates by </> 4 field theory calculations at fixed dimensions [14] 
and numerical Monte Carlo simulations [19] (see Table IV). 

The 2-d XY model (N = 2) is expected to follow the pattern of a Kosterlitz-Thouless 
critical phenomenon, whose critical region is characterized by a correlation length diverging 
exponentially with respect to r = (3 C — (3: £ ~ exp(6/r !T ) with a = 1/2. For this model the 
values of /,/ m (/3) and fi/ m (E) respectively at (3 C = 0.559(3) [29,30] and E c = 0.722(3) [29] 
are reported in Table II. Fig. 10 shows various non-defective /z/ m (/3) at N = 2, comparing 
them with the Monte Carlo results obtained recently by Kim [20] for correlation lengths: 
5 ^ £ ^ 70. The agreement is very good especially for PA's obtained using all available 
13 terms of the strong coupling series. PA's of the series in E turn out to be more stable, 
as shown in Table II and in Fig. 11, giving a perfectly consistent result for /*. Our final 
estimate is /* = 18.2(2) which is slightly larger than the Monte Carlo result /* = 17.7(2) [20] 
(this number has been obtained by taking only data for £ > 25 of Ref. [20] and taking into 
account the extra factor N in our definition (58)), but definitely consistent. 

When N > 3 the critical point moves to infinity making the determination of /* from 
strong coupling harder. For such models our analysis should be considered just exploratory 
due to the shortness of the available series, but as we will see the results look promising. In 
order to give an idea of the stability of our resummation procedure in this case, in Table III 
we report the values taken by fi/ m ((3) and fi/ m (E) for couples of (3 and E corresponding to 
a correlation length £ ~ 10. In Fig. 12 various /z/ m (/3) at N = 3 are drawn and compared 
with the Monte Carlo results of Ref . [20], obtained for correlation lengths 10 ^ £ ^ 120. The 
curves corresponding to different /z/ m (/3) are very close up to (3 ~ 0.5 (see also Table III). 
At (3 > 0.5 we observe that curves from different PA's become more and more stable with 
increasing n = I + m + 1, improving the agreement with the Monte Carlo data. Anyway, 
the agreement is quite good, even for (3 ~ 0.6, corresponding to £ ~ 100. Fig. 13 shows 
some approximants fi/ m (E) computed from the series in the energy. At E = 1 they give 
consistent results within 5-10%. Similar results are observed for N = 4, as shown in Fig. 14, 
where some /z/ m (/3) are plotted. Notice that at (3 = 0.6 £ ~ 25 [31]. 

In order to get an estimate of /* for both N = 3 and N = 4 we considered the values of 
fi/m((3) and fi/ m (E) at the largest values of (3 and E where they are still stable, i.e. (3 ~ 0.5 
and E ~ 0.6, which correspond to an acceptably large correlation length £ ~ 10. So, from 
data in Table III, our final estimate at N = 3 is g* ~ 19.8, with an uncertanty of few per 
cent, which compares very well with the Monte Carlo result /* = 19.6(2) (obtained by fitting 
all data reported in Ref. [20] to a constant) and 0(1/ N) calculation /* ~ 20.09. At N = 4 
our estimate is /* ~ 21.2 against /* ~ 21.35 coming from the 0(1/ N) calculation. 

Fig. 15 summarizes our 2-d results: it shows our strong coupling estimates of /* versus 
iV, comparing them with the available estimates of /* from alternative approaches: </> 4 
field theory, Monte Carlo and 1/N expansion techniques. There is a general agreement, in 
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particular the 0(1/ N) calculation /* ~ 8?r (l - ° 60 ^ 033 ) fits very well data down to N = 3. 
Furthermore we observe the linear approach of /* toward zero for iV — > — 2, similarly to the 
d = 1 case. 

Let's now consider 3-d 0(N) a models, which present a critical behavior at a finite (3 for 
all values af N. In order to check possible systematic errors we analyzed both the strong 
coupling series of f((3) and f(E). Table V shows a summary of the estimates of /* from 
the values of the approximants fi/ m (/3) at (3 C and fi/ m (E) at E c . Final estimates of /* from 
the analysis of the strong coupling series in (3 and in E are reported in Table VI. We recall 
that the errors diplayed in the Table VI are related to the spread of the PA results, while an 
estimate of the true systematic errors could only come from the comparison of results from 
different series associated to the same quantity. 

Fig. 16 shows typical curves of f(f3) obtained by [6/6] PA's (sometimes we used [7/5] 
PA's when the [6/6] ones were defective). The error bars displayed at (3 / (3 C = 1 show the 
spread of the estimates of /* from different non-defective [l/m] PA's. 

At large iV, typically N > 3, both series of f(f3) and f(E) give consistent results, which 
should be an indication of small systematic errors. As further check of our resummation 
procedure in the large- N region, we repeated our analysis at N = oo. We found that most of 
the approximants fi/ m ((3) constructed from the 13 th order series of A s ((3), if plotted in Fig. 6, 
would not be distinguishable from the exact curve. The analysis of the N = oo 13 th order 
series of A 3 (f3) and B 3 (E) would have given respectively /* = 50.25(6) and /* = 50.27(6) 
against the exact value /* = 16-7T = 50.2654.... Therefore everything seems to work fine 
at large N. On the contrary, at small N there are discrepancies between the analysis in (3 
and in E, which are definitely larger than the typical spread of the PA estimates of /* from 
each series. Such differences give somehow an idea of the size of the systematic errors of our 
analysis when applied to these values of N . 

In Table VI for comparison we give also the results from </> 4 field theory and 1/N expan- 
sion. Fig. 17 summarizes all available results for /*. There as a strong coupling estimate of 
/* we show the average of the results from the series f(f3) and f(E), while their difference 
is used as an estimate of the systematic error. 

At large iV, N > 8, there is a substantial general agreement: estimates from the strong 
coupling approach, 0(1/ N) calculation /* ~ 16-7T (l — 1 5 ^ r 601 ) , and </> 4 field theory differ at 
most by 1% to each other. At small iV, N = 0,1,2, our strong coupling estimates show 
relevant discrepancies with the field theoretical calculations, which are of the size of the 
differences between the results coming from the analysis of f((3) and f(E), and therefore 
they should be caused by systematic errors in the strong coupling analysis employed. Anyway 
such discrepancies are not dramatic, indeed they are at most 5% and decrease with increasing 
N. 

In conclusion we have seen that, in two and three dimensions 13 terms of the strong 
coupling series of Ad((3) and Bd(E) are already sufficient to give quite stable results, which 
compare very well with calculations from other techniques, such as </> 4 field theory at fixed 
dimensions, Monte Carlo simulations and 1/N expansion. Of course an extension of the 
series of / would be welcome, especially for two reasons: 

(i) to further stabilize the PA's in the asymptotically free models and obtain reliable 
estimates at values of (3 corresponding to large correlation lengths £ > 100, and moreover 
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check if the change of variable (3 — >■ E allows one to get a reliable strong coupling estimate 
of /* in the continuum limit E — > 1; 

(ii) to see if the apparent discrepancies at small N in 3-d with the more precise </> 4 field 
theory calculations get reduced. 

An extension of the series of f(f3) may also allow more accurate and flexible analysis, 
like differential approximants, which in general require many terms of the series in order to 
give stable results, and which could provide a better reconstruction of f(f3) from its strong 
coupling series, taking properly into account the confluent singularities, which should be the 
major source of systematic error in 3-d models at small N. 

VI. CONCLUSIONS 

We computed the dependence of the renormalized four-point coupling g r from the renor- 
malized mass m r and the bare coupling to 0(1 /N) for (9(iV)-invariant (</> 2 )^ theories (d < 4) 
in the symmetric phase. As a consequence we obtained expressions for the /3-function and 
its fixed point g* within the same approximation. 

We extracted an independent determination of g* from the strong coupling analysis 
of the 0(N) non-linear a models, which we performed for d = 2,3 and selected values 
of N in the whole range N > —2, applying resummation techniques both in the inverse 
temperature variable (3 and in the energy variable E. In two dimensions and for N sufficiently 
large (N > 3) in three dimensions we found a good agreement with the </> 4 fixed- dimension 
field theory estimates, and we could also check consistency with the 1/N prediction, thus 
seemingly indicating good convergence properties of the 1 /N expansion at least when applied 
to the above quantities. In three dimensions and for small iV, however, some discrepancy 
between resummations of the series in (3 and in E occurred, which we interpreted as an 
indication of systematic errors, and which was also reflected into a small disagreement with 
results presented in the literature and obtained with other techniques, like </> 4 field theory 
at fixed dimensions. Such discrepancy might be significantly reduced by knowing a few 
more terms in the strong coupling series, whose feasibility seems to be well within the range 
of present day strong coupling techniques. In our opinion improving the strong coupling 
analysis might lead to a determination of the fixed point value of the renormalized four-point 
coupling with a precision comparable to, or even better than, the best available results. We 
stress the crucial role played by the comparison of series in the variables (3 and E in order 
to estimate the relevance of systematic errors. 

ACKNOWLEDGMENTS 

It is a pleasure to thank A. J. Guttmann and A. D. Sokal for useful and stimulating 
discussions. 



21 



APPENDIX A 



By applying Eq. (35) to Eq. (24) and making explicit use of Eq. (36), we may obtain 
the following explicit representation of the 0(1/ N) contribution to the /3-function of 0(N) 
models in d dimensions: 
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where we have introduced the rescaled integration variable u = k/m r . By noticing that, 
according to its definition (37) 



> U(u) > 



(A2) 



it is easy to get convinced that all integrals appearing in Eq. (Al) are well-defined and 
finite as long as d < 4 and g r < g*. Moreover it is possible to perform a series-expansion 
of Eq. (Al) in the powers of g r , reproducing order by order standard perturbation theory 
results [17], and in particular to leading order: 



8$, 



(A3) 



for all values of d. 

For the sake of comparison in Figs. 18 and 19 we plot the function ^\g), where g has 
been defined as in Refs. [13,14,17] such that g* = 1 at N = oo (see also Sec. IV and Eq. (74)), 
respectively for d = 3 and d = 2. 



We recall once more that 
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and by definition we set 
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In Table VII we report all values in d = 1, 2, 3 such that (3 n > 10 3 . As a check of accuracy 
of the perturbative expansion we may employ the identity 
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6 + J2Pn = ~fi , (A6) 



n=l 



where fi was defined and evaluated in Sec. IV (cfr. Eqs. (65) and (72)). Notice that for d = 3 
the coefficients f3 n for n < 5 can also be extracted from the literature [17], and Eq. (A6) is 
already satisfied within 1% precision by the six-loop /3-function. We mention that, in the 
case d = 1, /3^\g) may actually be computed analytically, and the result is 

P (1 \g) = ^[g + ^r ( 648 - 3732 ^ + 5512 ^ 2 - 2m 9 3 - 330 ^ 4 - 27 9 5 ) 



Qg 



2 



( 3 + S) 



(81 - 6g + 1750g 2 - 1598# 3 + 509# 4 ) . (A7) 



The definitions (A4) and (A5) are obviously inappropriate in the limit d — > 4, in which 
case one may verify that 

— 8(3 f r - 9ftf r + 0(gt) (A8) 

where (3 — > j^z- Eq. (A8) in turn can be compared to the known perturbative evaluation 
around d = 4: 

m) = (d - A)g r + ^« - H3N N + 2 U) ^o9r + 0{gt) , (A9) 

finding complete agreement to 0(1 /N). 

It is conceivable to reinterpret the d — > 4 limit of Eq. (Al) in a nonperturbative sense 
by a principal-part prescription for the singularity occurring at the Landau pole u identified 
by the condition 

U(u) = . (A10) 

Work in this direction is in progress. 
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FIGURES 



FIG. 1. Feynman rules for the 1/N expansion. 
FIG. 2. Graphical definition of fundamental integrals. 
FIG. 3. Identities among Feynman graphs. 

FIG. 4. /* vs. N for 1-d 0(N) models. 

FIG. 5. vs. f3 for 2-d 0(oo) model. The dotted horizontal line represents the continuum 
value /* = 8x. 

FIG. 6. vs. f3 for 3-d 0(oo) model. The dotted horizontal line represents the continuum 

value /* = 16x. The dashed vertical line indicates the critical point f3 c = 0.252631.... 

FIG. 7. A/3 C = /3 c (oo) — )3 C (N) versus 1/N in 3-d models, as estimated by a strong coupling 
analysis. The dashed line represents the 0(1 /N) calculation (cfr. Eq. (67)). 

FIG. 8. AE C = E c (oo) — E C (N) versus 1/N in 3-d models, as estimated by a strong coupling 
analysis. The dashed line represents the 0(1 /N) calculation (cfr. Eq. (69)). 

FIG. 9. /6/6(/3) vs. PIPc f° r various values of N < 2 in two-dimensional models. 

FIG. 10. Some /j /m (/3) are plotted versus /3 for the 2-d XY model (N = 2). The vertical dotted 
lines indicate the critical point: f3 c = 0.559(3) [29,30]. Monte Carlo data from Ref. [20] are also 
shown. 

FIG. 11. Some fi/ m (E) are plotted versus E for the 2-d XY model (N = 2). The vertical 
dotted lines indicate the value of the energy at the critical point: E c = 0.722(3), estimated from 
Monte Carlo data [29]. Monte Carlo data from Ref. [20] are also shown. 

FIG. 12. Some /;/ TO (/3) are plotted versus f3 for the 2-d 0(3) model. Monte Carlo data from 
Ref. [20] are also shown. 

FIG. 13. Some fi/ m (E) are plotted versus E for the 2-d 0(3) model. Monte Carlo data from 
Ref. [20] are also shown. 
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FIG. 14. Some fi/ m (f3) are plotted versus f3 for the 2-d 0(4) model. 



FIG. 15. For 2-d models we plot /* vs. N obtained from our strong coupling analysis. For 
comparison field theoretical and Monte Carlo estimates are also shown. The dashed line represents 
the 0(1/N) calculation of g*. 

FIG. 16. vs. /3//3 c for various values of N in three dimensional models as obtained by a 

[6/6] PA (or [7/5] when the [6/6] one was defective). Error bars at /3//3 c = 1 show the spread in 
the determination of /* from all PA's considered. 

FIG. 17. /* vs. N from our strong coupling analysis in 3-d. For comparison field theoretical 
estimates are also shown. The dashed line represents the 0(1/N) calculation of /*. The dotted 
line indicates the value of /* at N = oo. 

FIG. 18. We plot (3^(g) vs. g for the three-dimensional case. 

FIG. 19. We plot (3^(g) vs. g for the two-dimensional case. 
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TABLES 



TABLE I. For some 2-d O(N) a models with N < 2 we report: the critical point, f3 c ; the 
estimate of the singularity of the PA's closest to the origin, /3 , which corresponds to a regular zero 
of A 2 (/3); the values of the approximants /;/ TO (/3) at the critical point. Asterisks mark defective 
PA's. 



N 


/3c 


A> 


5/5 


4/6 


5/6 


4/7 


6/6 


5/7 


4/8 


-1 


0.3145(1) 


-0.1315 


5.27 


5.32 


5.27 


5.28 


5.27 


5.31 


* 


-1/2 


0.3492(1) 


-0.1506 


7.87 


* 


* 


8.04 


7.85 


8.10 


8.00 





0.379052(1) [33] 


-0.1653 


10.51 


* 


* 


10.43 


10.54 


10.48 


* 


1/2 


0.408545(8) 


-0.1774 


12.60 


12.72 


12.68 


* 


12.66 


12.61 


12.63 


1 


0.4406867... 


-0.1878 


14.65 


14.72 


14.70 


14.67 


14.69 


* 


14.57 


3/2 


0.4804(1) 


-0.1969 


16.76 


16.76 


16.76 


16.60 


16.83 


16.83 


16.47 



TABLE II. For the 2-d XY model (N = 2) we give some details on the analysis of the series 
of /(/3) (first line) and f(E) (second line). We report: the critical point, f3 c (E c ); the estimate of 
the regular zero of A 2 (/3) (B 2 (E)) closest to the origin, f3 (E ); the values of the approximants 
fi/m{P) (fi/m(E)) at f3 c (E c ). Asterisks mark defective PA's. 



N 








5/5 


4/6 


5/6 


4/7 


6/6 


5/7 


4/8 


2 


/3c 


=0.559(3) [29,30] 


/3 =-0.2 49 


19.27 


19.44 


* 


18.71 


* 


* 


18.24 




E c 


=0.722(3) [29] 


E =-0.2179 


18.28 


18.46 


18.30 


18.35 


18.29 


* 


18.17 



TABLE III. We give some details of the strong coupling analysis of the series /(/3) (first line) 
and f(E) (second line) for two asymptotically free models: N = 3,4. We report: the estimate of 
the regular zero of A 2 (/3) (B 2 (E)) closest to the origin, f3 (E ); the values of the approximants 
fi/m((3) (fi/ m (E)) at a value f3 (E) corresponding to a correlation length £ ~ 10. The values of E 
and £ are taken from Ref. [32] for N = 3, and Ref. [31] for N = 4. Asterisks mark defective PA's, 



i.e. 


PA's 


having 


singularities 


for p < p (E 


<E). 














N 










5/5 


4/6 


5/6 


4/7 


6/6 


5/7 


4/8 


3 


/3o 


=-0.2188 


/3 =0.5 


11.05(1) 


20.3 


20.6 


* 


20.0 


* 


19.5 


19.8 




Eq 


=-0.2330 


E =0.60157 


11.05(1) 


19.9 


20.0 


19.9 


19.9 


19.9 


19.9 


19.8 


4 


/3o 


=-0.2305 


/3 =0.525 


10.32(3) 


21.8 


22.4 


* 


21.3 


* 


20.6 


21.0 






=-0.2456 


E =0.60089 


10.32(3) 


21.2 


21.3 


21.3 


21.3 


21.2 


21.4 


21.1 
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TABLE IV. For2-dO(iV) 


a models we report: 


the critical point fl c ; f* from our strong coupling 


analysis, /* c ; the 0(1 /N) calculation of /*, f£i N ; f* from </> 4 field theory at fixed dimensions, f£; 


f* from Monte Carlo simulations, /^ c . 




N f3 c 


/* 

J SC 


/l/jV /ft /mc 


-1 0.3145(1) 


5.29(3) 




-1/2 0.3492(1) 


8.0(1) 




0.379052(1) [33] 


10.51(5) 




1/2 0.40oo4o^oJ 


1*1 i2 O / CT \ 
12.00^0 J 




1 0.4406867... 


14.63(7) 


15.5(8) [14] 14.3(1.0) [19] 


3/2 0.4804(1) 


16.7(2) 




2 0.559(3) [29,30] 


18.2(2) 


17.7(2) [20] 


3 oo 


19.8(4) 


20.09 19.6(2) [20] 


4 oo 


21.2(5) 


21.35 


oo oo 




25.1327... 



TABLE V. For 3-d O(N) a models we present a summary of the analysis of the strong coupling 
series of /(/3) (first line) and f(E) (second line) at some values of N. We report: the critical point, 
/3 C (E c ); the estimate of the regular zero of A 3 (/3) (B 3 (E)) closest to the origin, f3 (E ); the values 
of the approximants /;/ TO (/3) (fi/ m (E)) at f3 c (E c ). Asterisks mark defective PA's. Errors due to the 
uncertainty of f3 c and E c are at most of the order of one in the last digit of the numbers reported 
(except for some cases where they are given explicitly). We mention that at N = 1 and N = 2 our 



strong coupling analysis led to E c = 0.332(3) for both. 



N 


/3c, E c 


/3o, Eo 


5/5 


4/6 


5/6 


4/7 


6/6 


5/7 


4/8 


-1 


0.19840(3) 


-0.109 


* 


11.2 


11.3 


10.9 


10.6 


10.7 


* 




0.350(5) 


-0.117 


9.3 


9.4 


9.7 


* 


10.1 


10.5 


10.3 





0.21350(1) [33] 


-0.134 


19.7 


19.8 


19.6 


19.8 


* 


19.4 


* 




0.333(5) 


-0.146 


* 


18.1(3) 


* 


18.4(2) 


18.3(2) 


* 


* 


1 


0.221652(4) [34] 


-0.149 


25.4 


26.0 


25.3 


25.8 


25.4 


24.8 


* 




0.3301(1) [36] 


-0.166 


24.3 


24.3 


24.3 


24.4 


24.4 


24.4 


24.4 


2 


0.22710(1) [35] 


-0.160 


29.4 


29.6 


29.4 


29.6 


29.4 


* 


29.6 




0.3297(2) [37] 


-0.180 


28.9 


28.9 


28.9 


29.0 


28.9 


28.9 


28.9 


3 


0.231012(12) [38,39] 


-0.168 


32.5 


32.3 


* 


32.4 


32.5 


32.4 


32.4 




0.331(3) 


-0.191 


32.2 


32.2 


32.2 


32.3 


32.3 


32.3 


32.3 


4 


0.2339(1) 


-0.175 


34.9 


34.5 


35.3 


34.6 


34.9 


34.7 


34.7 




0.333(2) 


-0.200 


34.8 


34.8 


34.7 


34.9 


34.8 


34.9 


34.8 


8 


0.2407(1) 


-0.19 


40.5 


39.9 


41.1 


40.2 


40.6 


40.4 


40.3 




0.334(1) 


-0.224 


40.5 


40.5 


* 


40.8 


40.6 


40.7 


40.7 


16 


0.2458(1) 


-0.20 


44.8 


44.3 


45.2 


44.6 


44.9 


44.7 


44.7 




0.3370(5) 


-0.246 


44.8 


44.8 


* 


45.1 


44.9 


45.0 


45.0 


24 


0.2479(1) 


-0.21 


46.5 


46.1 


46.8 


46.4 


46.6 


46.4 


46.4 




0.3379(3) 


-0.260 


46.4 


46.5 


47.0 


46.7 


46.6 


46.7 


46.6 


32 


0.2492(2) 


-0.22 


47.5 


47.1 


47.6 


47.3 


47.5 


47.4 


47.4 




0.3384(3) 


-0.268 


47.3 


47.4 


47.7 


47.5 


47.5 


47.5 


47.5 


48 


0.2502(1) 


-0.22 


48.4 


48.1 


48.5 


48.3 


48.4 


48.3 


48.3 




0.3390(3) 


-0.280 


48.2 


48.3 


48.5 


48.5 


48.4 


48.4 


48.4 


oo 


0.252731... 


-0.25 


50.24 


50.12 


50.31 


50.21 


50.27 


50.25 


50.24 




0.340537... 


-0.34 


50.15 


50.23 


50.36 


50.32 


50.26 


50.29 


50.26 
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TABLE VI. For 3-d 0(N) a models we report: /* as estimated by the analysis of the strong 
coupling expansion of f(P), / s * c/3 ; /* from the analysis of the series of f(E), /* cE ; the 0(1/N) 
calculation of /*, f£/ N ; f* from </> 4 field theory at fixed dimensions, f£ t . In Ref. [17] data of /* were 
reported without errors, and differences with Refs. [13,14] should be due to a different resummation 
procedure. 



JV 


Ac, ,3 


Ac,E 


A/JV 


/ft 


-1 


10.7(4) 


10.3(6) 









19.4(3) 


18.3(3) 




17.86(5) [14,13] 17.62 [17] 


1 


25.1(5) 


24.4(1) 




23.72(8) [14,13] 23.47 [17] 


2 


29.5(1) 


28.9(1) 




28.27(8) [14,13] 28.03 [17] 


3 


32.4(1) 


32.3(1) 




31.78(9) [14,13] 31.60 [17] 


4 


34.8(3) 


34.8(1) 


30.84 


34.41 [17] 


8 


40.4(4) 


40.7(1) 


40.55 


40.93 [17] 


16 


44.8(3) 


45.0(1) 


45.41 


45.50 [17] 


24 


46.5(2) 


46.6(2) 


47.03 


47.13 [17] 


32 


47.4(2) 


47.5(1) 


47.84 


47.94 [17] 


48 


48.3(1) 


48.4(1) 


48.65 




oo 


50.25(6) 


50.27(6) 


50.2654... 





TABLE VII. We report the values of /3 n , defined in Eq. (A5), in d = 1, 2, 3 such that /3 n > 10" 3 . 
Notice that in 2-d: Pi = f + ^ 2 - f ^'(1/3) = -10.33501055. 





d= 1 


d= 2 


d = 3 


Pi 


388 
27 


-10.33501055 


164 
27 


Pi 


1187 
108 


5.00027593 


1.34894276 


Ps 


335 
162 


-0.08884297 


0.15564589 


Pa 


10001 
46656 


-0.00407962 


0.05123618 


Ps 


605 
11664 


0.00506747 


0.02342417 


Pe 


20045 
1119744 


0.00491122 


0.01264064 


Pr 


38671 
5038848 


0.00377364 


0.00757889 


Ps 


1231807 
322486272 


0.00281096 


0.00489401 


Ps 


21367 
10077696 


0.00211235 


0.00334024 


PlO 


89062753 
69657034752 


0.00161697 


0.00237987 


Pll 


28651973 
34828517376 


0.00126267 


0.00175481 


Pl2 




0.00100476 


0.00133070 


Pl3 




0.00081329 


0.00103290 


PlA 






0.00081770 
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